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We discuss the (non)-restoration of global and local symmetries at high tempera- 
ture. First, we analyze a two-scalar model with Z2 x Z2 symmetry using the exact 
renormalization group. We conclude that inverse symmetry breaking is possible in 
this kind of models within the perturbative regime. Regarding local symmetries, 
we consider the SU (2) (g) U (1) gauge symmetry and focus on the case of a strongly 
interacting scalar sector. Employing a model-independent chiral Lagrangian we 
find indications of symmetry restoration. 
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1 Review of Symmetry Non-restoration and Inverse Symmetry 
Breaking 

It has been known for a long timeB, that simple multi-scalar models can ex- 
hibit an anti-intuitive behaviour associated with more broken symmetry as 
the temperature is increased. We refer to this behaviour as inverse symmetry 
breaking or symmetry non-restoration, depending on whether the symmetry 
is exact or not at zero temperature. This phenomenon may have remarkable 
consequences for cosmology □. 

In renormalizable supersymmetric theories internal symmetries are always 
restored at high temperature, and there are arguments of why this should 
be also the case even when non-renormalizable terms are considered; symme- 
try non-restoration cpuld be possible though in the presence of non-vanishing 
background chargesu. 

Doubts have been raised on the validity of these results, since they are 
based on the one-loop approximation to the finite temperature effective po- 
tential, which is known to be unreliable for the discussion of many aspects of 
phase transitions. Different techniques are being actively applied to improve 
the one-loop approximation, leading quite often to contradictory results. The 
effect of next-to-leading-order contributions within perturbation theory has 
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been investigated^, through the study of gap equations, which are equivalent 
to a resummation of the super-daisy diagrams of the perturbative series. Large 
subleading corrections have been identified, which lead to a sizeable reduction 
of the parameter space where inverse symmetry breaking occurs. The question 
of inverse symmetry breaking has also been studied through the use of the 
renormalization group and a variational approach 0, with similar conclu- 
sions. Contrary to the results of the above studies, a large- N analysis seems 
to indicate that symmetry is always restored at high temperature However, 
the validity of this claim has recently been questioned A finite-lattice cal- 

Slation also supports symmetry restoration at sufficiently high temperature 
. Although the relevance of this result for the continuum limit is not clear, 
| a_Monte Carlo simulation in 2+1 dimensions seems to support this conclusion 



2 Inverse Symmetry Breaking and the Renormalization Group 

In this section we consider the simplest model that exhibits inverse symmetry 
breaking: a two-scalar model with Z2 x Z2 symmetry. The tree-level potential 
is given by 

Vtr(<h.><h) = 2 m l^l + ^™2^2 + jAl^ + ~A 2 </>| ~ ^Al2^0| • (1) 

This potential is bounded for Ai, 2 > and 

AiA 2 > \\ 2 . (2) 

In the high temperature limit, |m»| <C T, Xhe thermal correction to the above 
potential at the one-loop level is given byl!3 

Ay T (</» 1 ,</» 2 )~^-[(3A 1 -A 12 )0? + (3A 2 -A 1 2)0l] + •■■ (3) 
For the parameter range 

3Ai - A12 < , (4) 

which can be consistent with the stability condition of eq. (||), the thermal 
correction for the mass term of the <f>i field is negative. Notice that the stability 
condition (0) does not allow both mass terms to be negative. If the system 
is in the symmetric phase at zero temperature with m\ 2 > 0, there will be a 
critical temperature T c 2 r = 12mf/(Ai2 — 3Ai) above which the symmetry will 
be broken. If the system is in the broken phase at T = 0, the symmetry will 
never be restored by thermal corrections. 
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Our aim is to discuss the above scenario in the context of the Wilson 
approach to the renormalization group. The main ingredient in this approach 
is an exact flow equation that describes how the effective action of the system 
evolves as the ultraviolet cutoff is lowered. We consider the lowest order in a 
derivative expansion of the effective action, which contains a general effective 
potential and a standard kinetic term. At non-zero temperature this approach 

an be formulated either in the imaginary-timeEj or in the real-time formalism 
. In the latter formulation, the evolution of the potential lowering the cutoff 
scale A is given by the partial differential equationtil 



-T^Tr^l,,, 



1 — exp 



(5) 



where 



[M 



a 2 y A (0i,^) 



,J = 1,2 



(6) 



The initial condition for the above equation, at a scale Ao 3> T, is the renor- 
malized effective potential at zero temperature. We consider small quartic 
couplings, so that the logarithmic corrections of the zero-temperature the- 
ory can be safely neglected. The initial condition for the evolution is a zero- 
temperature potential given by eq. ([!]) . Integrating the evolution equation (0) , 
we obtain the non-zero-temperature effective potential in the limit A — > 0. 

Finding the solution of eq. (|5|) is a difficult task. An approximate solution 
can be obtained by expanding the potential in a power series in the fields. In 
this way the partial differential equation (^) is transformed into an infinite 
system of ordinary differential equations for the coefficients of the expansion. 
This system can be solved approximately by truncation at a finite number of 
equations. That is, the potential is approximated by a finite-order polynomial. 
As a first step, we follow this procedure and define the running masses and 
couplings at the origin 



m? 2 (A) 



d 2 v A 



,Ai,a(A) 



01,2=0 



ld 4 V A 



G i 



,Ai 2 (A) 



01,2=0 



1 d 4 V A 
2 



01,2=0 

(7) 

The corresponding evolution equations can be obtained by differentiating eq. 
(p|) and neglecting the higher derivatives of the potential. We find 



d_ 
*BA' 



A^—m 2 



A 9 , 

dA ' 2 



1,2 — ~ 6Ci ,2-^1,2 + 2C24A12 

— lS-Di^A 2 2 — 2Z?2,iA^2 



(8) 
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with 



a|-A 12 = -6D1A1A12 - 6 D 2 X 2 \ 12 + 8 °\ % \\ 2 , 
oA m\ z — mi 



A 3 N(m h2 ) dC l>2 2 2 2 

C l>2 = — — , D l,2 = a _ 2 ■ ^l,2= A + m l,2> ( 9 ) 



47T CJx 2 



'1,2 



and JV(a;) = [exp(a;/T) — 1] the Bose-Einstein distribution function. For 
^1.2 <C 7 1 we have 

A 3 T 

C\a — > — r t „ — o , (10) 
An 2 A 2 + mf 2 

and the above equations agree with those considered in ref. i in the same limit. 
For 071,2 3> T there is no running, because of the exponential suppression in 
the Bose-Einstein function. 

We have solved numerically the system of equations (^J) and determined the 
range of zero-temperature parameters that lead to inverse symmetry breaking. 
In fig. 1 we present the results for a zero-temperature theory with positive mass 
terms m 2 (Ao) = to 2 (Ao) and A2(Ao) = 0.3. The temperature has been chosen 
much higher than the critical one (T = 500mi(Ao)). The system (||) has been 
integrated from Ao T down to A = 0, where the thermally corrected masses 
and couplings at non-zero temperature have been obtained. A negative value 
for the mass term m\ at A = has been considered as the signal of inverse 
symmetry breaking. This has been achieved in the region above line (a) in fig. 1. 
We also plot the stability bound of eq. @ (the allowed range is below line (b)), 
and the perturbative prediction for the range that leads to inverse symmetry 
breaking (above line (c)). The phenomenon of inverse symmetry breaking 
is confirmed by our study, in agreement with ref. Ej, where the imaginary- 
time formulation of the renormalization-group approach has been used. We 
observe that the renormalization-group treatment eliminates a large part of 
the parameter space allowed by perturbativa-theory, in agreement with the 
results obtained by solving the gap-equations □. 

The reliability of our conclusions crucially depends on whether the solution 
of the system of truncated equations (0) provides an approximate solution to 
the full partial differential equation (H). We have checked that tcz^numerical 
integration of eq. (||) through the algorithms discussed in ref. Due to 
limitations in computer time, we restrict our discussion of eq. along the 
4>i axis, which is the direction of expected symmetry breaking for our choice 
of couplings. We approximate the potential by the expression 

V A (&,&) = Wi) + ~m|(A)0| + ^X 2 (A)4 - ^\ 12 (K)4>\4>1 . (11) 
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Figure 1: Parameter space that leads to inverse symmetry breaking (A2 = 0.3). 



The evolution of m| (A) , A2 (A) and A12 (A) is determined through the truncated 
eqs. (||). However, the full <f>\ dependence is preserved through the numerical 
integration of eq. (||), with the eigenvalues of the mass matrix fA\ given by 

Wi]^^^ ^d [Mi] 2 = mI(A)-A 12 (A)^. (12) 

This treatment permits a reliable study of the order of the symmetry-breaking 
phase transition; we have found that it is governed by the Wilson-Fisher fixed 
point of the onepficalar three-dimensional theory, resulting in a second order 
phase transition Q. 

3 SU(2) (g) U(l) gauge symmetry with strongly coupled Higgs sector 

A natural question to ask is whether the standard model gauge symmetry 
SU(2) ® U (1) could remain broken at high temperatures. It is well known that 
the symmetry is restored in the minimal standard model, so we consider its 
simplest extension with two Higgs doublets. We find that due to the positive 
and large contribution of the gauge bosons to the scalar thermal masses, it 
is not possible to attain a negative mass term (needed for symmetry non- 
restoration) within the perturbative range of the scalar couplings. 

A strongly coupled Higgs sector implies (at least naively) heavy physical 
scalar particles, which can be effectively removed from the physical low-energy 
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spectrum. To study the behaviour of the gauge symmetry in this case we use 
an effective Lagrangian which keeps only the light degrees of freedom, namely 
the gauge and Goldstone bosons together with the fermions. The resulting 
chiral Lagrangian is a non-renormalizable non-Jinear sigma model coupled in 
a gauge invariant way to the Yang-Mills theoryES. At lowest order, it is model 
independent. 

We propose to use the gauge boson magnetic mass as an indicator of sym- 
metry (non)-restoration. ft is defined as the transverse part of the correspond- 
ing self-energy, HT(0,fc), on-shell. A perturbative computation shows that it 
is exactly equal to zero at one loop in an unbroken gauge theory. For unbro- 
ken non-Abelian gauge theories, such as QCD, higher orders in perturbation 
theory suffer from infrared divergences, and a magnetic mass of order g 2 T is 
expected to be generated non-perturbatively. In spontaneously broken gauge 
theories, such as the standard electroweak model and its extensions, no such 
divergences are present. Thus, we expect that even in perturbation theory the 
magnetic mass (as computed in the broken phase) will show a tendency to 
vanish at high enough temperatures, whenever symmetry restoration occurs. 

We have calculated the thermal gauge boson self-energies at one loop and 
leading order, 0(T 2 ), from which we obtain the magnetic masses 



M^ mag = 5 2 ^, (13) 
v(T) 2 



Mi mag = (g A +g' 2 )^, (14) 

with 



v(T) 2 = v 



T 2 

6v 9 



where the fact that the gauge couplings are not renormalized at this order has 
been used. We conclude thus that in models with strongly interacting Higgs 
sector the spontaneously broken SU(2) ® U(l) gauge symmetry tends to be 
restored when the system is heated. Notice though that our calculation is only 
valid for temperatures below the electroweak scale. Details of the calculation 
and the approximations involved can be found inEZI. 

It is worth to remark that the thermal corrections to v coincide with those 
of the pion decay constant F„ in the non-linear sigma model E3. That is, at 
one loop and leading order T 2 all the temperature corrections to the scalar 
v.e.v. are due to the would-be Goldstone bosons; fermions and gauge boson 
transverse degrees of freedom will only contribute at higher order. 
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